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Abstract A-122 pnictides {A an alkaline earth ele¬ 
ment) possess striped antiferromagnetic states of or¬ 
thorhombic (OR) symmetry and nonmagnetic states of 
tetragonal (T) symmetry. Based on total energy cal¬ 
culations, we show that the T = OK pressure-driven 
transition from OR to T states occurs at a pressure, 
Pff, where the tetragonal enthalpy drops below the or¬ 
thorhombic. The OR state remains metastable up to a 
higher pressure. Pm > Ph- We find anomalies in mag¬ 
netism and orthorhombicity as P —> Pm, and a trend 
towards elastic instability. 

Keywords 122-Pnictides • Elastic instability • 
Metastable OR • Pressure 

1 Introduction 

The 122 pnictides, AFe 2 As 2 [A = Ca, Sr, Ba), dis¬ 
play structural, magnetic, and superconducting phase 
transitions upon doping or applied pressure [TJ[11[31I1]. 
At ambient pressure, Ca-, Sr-, and Ba-122 compounds 
transition from the high-temperature tetragonal (T) 
phase to a low temperature orthorhombic (OR) phase, 
striped along the a-axis and antiferromagnetically (AFM) 
ordered along the c-axis, at Neel temperature Tjq ^ 
170K, 205 K, 140 K respectively [hllHllT] . The T-OR 
transition temperature decreases with applied pres¬ 
sure. 
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At low-T, under hydrostatic pressure, the 122-pnictides 
exhibit rich universal behavior, albeit with interesting 
differences. While Ca-122 transitions from the OR-AFM 
phase to a collapsed tetragonal (cT) phase with strongly 
decreased c-axis at P ~ 0.35 GPa [S], Sr-122 and Ba- 
122 transition to a T phase with only slightly decreased 
c-axis, at P ^ 4.4 GPa and 10.5 GPa respec¬ 

tively mm- These eventually transition to a cT phase 
at higher pressures. 

Our previous work [T3], based on first principles 
density functional theory calculations, showed that the 
pressure-driven T = OK transition from OR-AFM to 
tetragonal (cT or T) can be understood as a l®*^-order 
enthalpic transition that occurs at a pressure we de¬ 
note Ph- Subsequent work [Hill], also based on first 
principles total energy calculations, proposed that the 
T - cT transition and anomalies in lattice parameters 
and elastic properties, observed experimentally at fi¬ 
nite temperatures, can be interpreted as arising from 
proximity to T = OK Lifshitz transitions, at pressures 
denoted as Pl, wherein pressure-driven changes in lat¬ 
tice constants alter the Fermi surface topology in these 
materials. Finally, our work revealed the existence of a 
metastable OR-AFM phase that remains mechanically 
stable until it transitions discontinuously to the tetrag¬ 
onal phase at a pressure we denoted Pm- 

Here we explore A-122 pnictides in further detail, fo¬ 
cusing on properties of the metastable OR-AFM phase. 
We study the pressure dependence of relative enthalpy, 
magnetism, and lattice parameters, which all vanish 
discontinuously at Pm- Additionally, elastic stability 
vanishes continuously as P —f Pm- 
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2 Methods 

Our calculations follow methods outlined previously [la 
E], using projector-augmented wave density functional 
theory in the PBE generalized gradient approximation |16] 
as implemented in VASP m Owing to the previously 
reported m sensitivity to calculational details, in the 
present work we take low Fermi surface smearing of 
tr = 0.05 eV, with a high density 15x15x7 fc-point grid 
for our 20-atom OR-AFM cells and matching 20-atom 
•\/2 X >72 tetragonal supercells. We cut off our plane- 
wave basis at 350 eV which is 30% greater than the de¬ 
fault for Fe-based alloys. Additional accuracy of forces 
is achieved though Fourier transform grids of full res¬ 
olution, additional support grid for the evaluation of 
augmentation charges, and application of projectors in 
reciprocal space. Self-consistency of electronic structure 
is achieved when subsequent iterations change the en¬ 
ergy by less than 5 x 10“^°eV. Relaxations are consid¬ 
ered converged when forces drop below 10“®eV/A and 
pressure matches its target to within 10“‘*GPa. Relax¬ 
ations are halted and restarted repeatedly to remove 
errors connected with volume-dependence of the basis 
sets. To calculate elastic constants we employ central 
finite differences using a small step size of 0.010 A, be¬ 
cause the stress-strain relation can be strongly nonlin¬ 
ear. 

In order to locate the end point of metastability. 
Pm-, we increment the pressure in steps from low to 
high, taking the relaxed structure at one pressure as 
the starting point for the next. Successively smaller 
steps are employed as the threshold pressure is ap¬ 
proached. Relaxation becomes intractable very close to 
Pm, because relaxation times diverge and error sensi¬ 
tivity grows. Hence we locate Pm only to within 0.1 
GPa. 


3 Results 

To distinguish stable OR-AFM, metastable OR-AFM 
and T states, we monitor the enthalpy difference, 

= AE + PAV = Tftetra ~ ^I^ortho) (1) 

the iron magnetic moment M, and orthorhombicity 
a — b 


We also monitor the elasticity tensor [18], which for 
orthorhombic symmetry is 
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Strains of Voigt type 1,2, and 3 correspond to stretch¬ 
ing, respectively, along the x, y and z axes, while type 
4, 5, and 6 correspond to shears of types yz, zx and xy. 
Enthalpy, magnetism and orthorhombicity, and elastic¬ 
ity are discussed in turn in the following subsections. 

3.1 Enthalpy 

As described in detail in our previous work[T^. and 
shown in the top panels (a) of Figs. 1-3, each of the 
three A-122 compounds undergoes a l®*-order enthalpic 
transition from the OR-AFM phase to a non-magnetic 
tetragonal phase (cT or T) where AH crosses 0 at pres¬ 
sure Ph- However, a thermodynamically metastable OR- 
AFM state persists up to a higher pressure, Pm > Ph- 
This metastable state has higher enthalpy than the 
competing nonmagnetic tetragonal state {i.e. AH < 0), 
but it retains mechanical stability. Thus it constitutes 
a stable fixed point of conjugate gradient relaxation, 
considered as a dynamical system whose attractors are 
local energy minima. Beyond Pm, this fixed point van¬ 
ishes. Relaxation carries the system towards an alter¬ 
nate fixed point, which turns out to be a nonmagnetic 
tetragonal state that is both thermodynamically and 
mechanically stable. Hence, AH vanishes discontinu- 
ously for P > Pm, as can be seen in panel (a) of each 
figure. 

Because the Lifshitz transition pressure Pl is nega¬ 
tive for Ga-122, while the metastable limit Pm is posi¬ 
tive, we find that Pl < Pm for Ca-122. Gonsequently, 
the metastable orthorhombic state relaxes to the col¬ 
lapsed tetragonal state with a low c-axis for P > Pm- 
In contrast, for Sr-I22 and Ba-I22, we find that Pm < 
Pl, so the metastable orthorhombic state relaxes to 
the uncollapsed tetragonal state with large c-axis for 
P > Pm, followed by a second transition (the Lifshitz 
transition [14]) to the collapsed state with low c-axis at 
higher pressures Pl > Pm- 

Values of Ph and Pm grow monotonically as we 
move down the alkaline earth column of the periodic 
table, from Ca through Sr to Ba. However, the frac¬ 
tional interval of metastability {Pm — Ph)/Ph dimin¬ 
ishes from a high of 700% for Ga-122 to a low of 5% for 
(2) Ba-122. 
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3.2 Magnetism and Orthorhombicity 

The tetragonal states (T and cT) lack magnetism and 
orthorhombicity. Thus, M and 6 ab vanish discontinu- 
ously at Pm as can be seen in panel (b) of each figure. 
They approach finite limiting values as P approaches 
Pm from below, while they vanish identically for all 
P > Pm- Enlarging the data near the discontinuity 
(not shown) reveals nonanalyticity in which M and 6 ab 
approach their limiting values with a singularity con¬ 
sistent with a square root power law. The similar be¬ 
havior of M and Sab reflects the coupling of magnetism 
and lattice parameters. 

3.3 Elasticity 

Notice that C is block diagonal, with a symmetric 3x3 
block for i,j < 3, and three diagonal entries for i = j > 
4. Thus nine elastic constants are required to describe 
the elastic response in the OR-AFM state. Tetragonal 
symmetry can be regarded as a limiting case of or¬ 
thorhombic symmetry for which orthorhombicity van¬ 
ishes {Sab = 0) and C obeys the identities Cii=C 22 , 
Cia = C 23 and C 44 = C 55 . Values of are plotted in 
panel (c) of each figure. 

Elastic stability requires C to be positive definite, 
that is, its eigenvalues {Xk} must be positive. This de¬ 
mands the matrix elements Cij obey a set of six con¬ 
straints, two of which are nonlinear [19] . As an alterna¬ 
tive test of positivity, we plot the eigenvalues in panel 
(d) of each figure. The eigenvalues are ordered so that 
Afc with fc = 1 — 3 denote the eigenvalues of the 3x3 
block (in increasing order), while /c = 4 — 6 simply re¬ 
peat the diagonal entries Ckk- When elastic stability is 
violated a crystal will spontaneously deform through a 
combination of strains characterized by the eigenvector 
of the vanishing or negative eigenvalue. 

All Cij exhibit anomalies as P —)■ Pm- The elastic 
constant C 33 in particular shows a strong dip for all A- 
122. Low C 33 is in danger of violating the Born criterion 
derived from positivity of the determinant of the 3x3 
block [T^ . 

^ . CnCia + C22C?3 + 2 C 12 C 13 C 23 

<-"33 > --7^2-’ 

U 11 O 22 — U 32 

which is equivalent to the requirement Ai >0. Near 
Pm we find that Ai « 0 to within our uncertainty 
that is limited by our difficulty in relaxation and our 
finite difference numerical differentiation. Ai tends to¬ 
wards 0 with approximately a square root singularity. 
The eigenvectors vi corresponding to eigenvalues Ai are 
superpositions of strains of types 1-3. The associated 
deformation stretches the c-axis while shrinking the a- 


and 6 -axes (or vice-versa). Orthorhombicity is evident 
in the differences between C13 and C23 which grow to 
around 50% as P ^ Pm, even while Cn « C 22 and 
C44 « C33 are nearly consistent with tetragonality. 


4 Discussion 

Relaxation of a crystal structure can be considered as 
a nonlinear dynamical system with dissipative dynam¬ 
ics [5D] . Coupling of the crystal atomic structure to the 
electronic band structure via the Kohn-Sham equations 
of density functional theory make the system effectively 
infinite dimensional. Such dynamical systems typically 
exhibit multiple fixed points, and these fixed points can 
bifurcate as parameters are varied. In the present case, 
we are confronted with the seeming disappearance of a 
stable hxed point corresponding to the OR-AFM state, 
as P ^ Pm, which would seem to violate local conser¬ 
vation of Poincare index . 

The mystery could be resolved if the stable fixed 
point were to annihilate with an unstable hxed point, 
such as occurs in a saddle-node bifurcation. The square- 
root singularities observed in magnetization, orthorhom¬ 
bicity, and elasticity are typical of such a bifurcation. 
A concrete demonstration of this hypothesis would re¬ 
quire observing the unstable hxed point. Unfortunately 
this is not readily achievable, as our conjugate gradient 
dynamics is optimized to seek stable hxed points. 

Considering that the numerical value of the thresh¬ 
old pressure Pm is sensitive to Fermi surface smear¬ 
ing tr, it is interesting to ask if the loss of metastabil¬ 
ity could be caused by a Lifshitz transition (change in 
Fermi surface topology). Indeed, we do hnd bands ap¬ 
proaching to within a of the Fermi energy as P ^ Pm 
for all three pnictides. However, the large unit cell with 
low symmetry and s, p, and d elements inevitably cre¬ 
ates many bands, leading to a proliferation of Lifshitz 
transitions, including several close to Pm, making it 
hard to assess the relevance. This raises the possibility 
that superconductivity which, in Sr- and Ba-122, oc¬ 
cur in the vicinity of Pm, may be correlated with these 
Lifshitz transitions [Hj. 

A pronounced dip in C 33 is indicative of soften¬ 
ing of associated phonon mode(s) as the system trans¬ 
forms from OR-AFM to a T phase. This is reminis¬ 
cent of Martensitic transitions in some of the A15 com¬ 
pounds |23], except that here this occurs as a function 
of pressure at T = 0. While the underlying mechanism 
for superconductivity in pnictides is believed to be mag¬ 
netic in origin, we conjecture that phonon mode soften¬ 
ing, and strong magneto-elastic coupling that our cal¬ 
culations suggest, may be consequential to the pressure 
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Fig. 1 Data for Ca-122. (a) Enthalpy difference ziiJ (Eq. |Q, 
units meV/atom), (b) magnetization M (units fin) and di¬ 
mensionless orthorhombicity Sat (Eq. ([^), (c) elastic con¬ 
stants Cij (units GPa), and (d) eigenvalues (units GPa) 
all as functions of pressure P (GPa). Plotting symbols indi¬ 
cate pressures where calculations were performed. Enthalpic 
transition at Ph and metastability limit Pm are marked. 



and temperature range over which superconductivity 
occurs. 
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